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We analyse the ground state of spinless fermions on a lat- 
tice in a weakly disordered potential, interacting via a nearest 
neighbour interaction, by applying the self-consistent Hartree- 
Fock approximation. We find that charge density modulations 
emerge progressively when r^ > 1, even away from half-filling, 
with only short-range density correlations. Classical geome- 
try dependent magic numbers can show up in the addition 
spectrum which are remarkably robust against quantum fluc- 
tuations and disorder averaging. 

PACS Numbers: 73.20Dx, 73.23Hk 

The interplay of disorder and interactions in two di- 
mensional Fermi systems is currently a central problem 
in condensed matter physics. Mesoscopic systems pro- 
vide a unique forum for analysing ground state proper- 
ties as it is possible to access the regime kT ^ A, where 
A is the mean single particle level spacing. Examples 
include the study of low temperature persistent currents 
and magnetic response in small quantum rings and dots, 
and low bias measurements of the d.c. response ^Q and 
capacitance [[7[ of weakly coupled quantum dots. The 
latter experiments made it possible to access directly the 
energy differences iin between ground states of N and 
A^ - 1 particles: fiN = E{N) - E{N - 1). The addition 
spectrum, ^^ (5(/i— /li), depends sensitively on the nature 
of the mesoscopic ground state [|| . 

Resonant tunnelling measurements of the addition 
spectrum [H-pf have resulted in some interesting observa- 
tions. Whilst the mean peak spacings are well described 
by the constant interaction model B, the fluctuations 
are not. In \^M the averaging was carried out over TV, 
whereas in IpETthe results were also averaged over sam- 
ple geometry. The experimental data 0] indicates the 
existence of atypical addition spectrum spacings at cer- 
tain values of TV, suggesting that averaging over disorder 
may not be equivalent to averaging over iV (the ergodic- 
ity principle |10| is violated). Theoretical and numerical 
studies [p|,0-d8J attempted to address various aspects of 
the problem. 

In the capacitance experiments of Ref. [uj, the mea- 
sured addition spectra display bunching, an indication 
that the Coulomb blockade becomes negative between 
one or more consecutive electron addition peaks; in the 
experiment, these peaks then coalesce. Such bunching is 
in direct conflict with the naive picture combining the 
constant interaction model with ergodic effective single- 
particle wavefunctions. It has been shown ||l^ that a clas- 
sical charge model can reproduce many of the observed 



effects, but there is currently no quantum mechanical ex- 
planation as the experiments are carried out at densities 
considered too high to form a Wigner solid (in the case 
of a Coulomb bare interaction see ref. ||2^], for a short- 
ranged potential one might expect a Wigner solid to be 
less stable). 

Motivated by these experiments, but not attempting 
to reproduce specific details thereof, we analyse the na- 
ture of the ground state by applying the self-consistent 
Hartree-Fock (SCHF) approximation. We are thus able 
to go beyond the random phase approximation (RPA) 
with perturbation theory (valid for large dimensionless 
conductance g, and small r^ [ [l2[[li|| ), whilst consider- 
ing larger systems than is feasible by exact methods. 
Experimental values of r^ > 1 have indeed been re- 
ported H. Starting from a noninteracting model, we 
find that as the interaction strength is increased such 
that Ts ^ 1 (but still too weak to form a Wigner solid), 
the electron gas crosses over to a regime where: (i) there 
are significant spatial density modulations; (ii) density- 
density correlation functions seem to saturate, defining 
only short range order; (iii) the addition spectrum be- 
comes strongly N dependent, with magic numbers for 
which A2(A^) = iiN+i ~ Mw exhibits sharp maxima that 
coincide with the related classical charge model. 

Recent reports of exact numerical studies ]2l|,p|, p2| , p2[ 
have emphasised that the properties of quantum dots 
which are not reproduced by effective single-particle 
random-matrix-like theories are associated with the 
emergence of short range correlations. Our results sup- 
port this claim, but the main thrust here is related to 
(iii): some deviations from RPA behaviour have a direct 
classical electrostatic counterpart. The signature of the 
latter is not totally washed out by quantum fiuctuations 
even far from the Wigner crystalisation threshold. 

We consider the following tight binding Hamiltonian 
for spinless fermions with periodic boundary conditions: 
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where (ij) denotes pairs of nearest neighbours, Wi is the 
random on-site energy in the range [—1^/2,1^/2], and t 
the hopping matrix element. All lengths are measured 
in units of the lattice constant a, so that U = e^/a and 
t = h'^/2ma^. For low filling (i.e. a parabolic band) 
we find rs = U/{ty/4^) H where v = N/A is the fill- 



ing factor, for N electrons in an area A. For the non- 
interacting system we find g = kpl/2 = 96TTi'{t/W)'^ by 
applying the Born approximation (valid for 1 ^ g ^ A) ; 
kp = \f^Kv and I is the elastic mean free path. In the 
capacitance measurements 0, the 2d quantum dots was 
sandwiched between a metallic source (heavily doped n+ 
GaAs) and drain (Cr/Au), separated, at distances com- 
parable with the mean particle separation, by tunnel bar- 
riers. To account for external sources of screening (taken 
as half planes), one can insert a bare interaction between 
electrons in the dot that is dipolar (1/r'^) at distances 
greater than the dot to gate separation when there is only 
one close gate, and in the case of two close gates with ex- 
ponentially small long range interactions |24|] . Here we 
model such interactions with a nearest-neighbour pair 
potential. 

The ground state is obtained in the SCHF approxi- 
mation, over a range of densities and disorder strengths 
at zero magnetic field. The generalised inverse partici- 
pation ratio (GIPR) is then calculated according to the 
following definition: 
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where p(\:i) denotes the expectation value of the total 
density at the lattice site i. The angle brackets corre- 
spond to an average over the disorder ensemble. The 
GIPR provides a convenient measure of the degree of den- 
sity modulation: in the limit of a perfectly flat density 
profile it takes the value unity, and increases for a mod- 
ulated density. The maximal value that can be obtained 
for the GIPR occurs when all the charge is concentrated 
on only N sites, in which case X = Xjv. 




FIG. 1. The GIPR is plotted for a range of disorder values 
as a function of the interaction strength U . The system is 
a 16 X 16 lattice with n.n. (solid) and Coulomb (dashed) 
interactions. We do not average over disorder, but remark 
that the GIPR is self averaging. Here v = 1/4, implying 
r, K. .^&U/t. 

The GIPR is plotted for a range of disorder strengths 
in Fig. 1. For v = 1/4 it increases rapidly between U k, 1 
and [/ « 5 depending on the disorder strength, then gives 



way to a weak interaction dependence for t/ > 5. For 
comparison we also plot results for an identical system 
but with bare Coulomb interactions, such that the in- 
teraction potentials are both equal to U between nearest 
neighbour sites: the relative rapidity of the increase of X 
for nearest neighbour interactions is clear. The increase 
in the GIPR signals an increase in the spatial modula- 
tion of the total electron density, we shall refer to the 
increased density modulation at finite [/ as a charge den- 



sity modulation (CDM) |25|. 

At zero interaction we find X — 1 ^ l/g for large g 
(not shown). Within our numerical accuracy we were 
unable to find a consistent size dependence in the GIPR, 
suggesting that disorder is the dominant mechanism con- 
trolling the small to large U cross-over, as seen in fig. 1 




FIG. 2. The density-density correlation function (see text) 
shows increasing short ranged correlations as U/t is increased. 
Here W — 4, the lattice is 16 * 16. Here v — 1/4, implying 
rs « .56U/t. The dotted line shows approximately one stan- 
dard deviation due to sampling, assuming that all correlations 
have vanished. Inset: at half filling we find no decay of the 
correlation function with distance for U > 1 

The GIPR yields no information on the spatial struc- 
ture of the ground state, for which we evaluate the 
density-density correlation function defined as 



C{r) = (p(r)p(0)),/(p(0)2 



(3) 



The subscript c indicates that only connected averages 
are included, and here, due to the homogeneity of the 
disorder averaged potential, the correlation function de- 
pends only on the vector separation r. We only consider r 
to be directed along a lattice vector (1,0). A typical result 
for the correlation function is plotted in Fig. 2. As the 
interaction strength is increased, short range correlation 
develop, and then saturate. The underlying square lat- 
tice excludes the possibility of observing incipient Wigner 
crystal fluctuations, which possess the symmetry of a tri- 
angular lattice. 

Comparing figures 1 and 2, one can see that the short 
ranged correlations develop over the same range of inter- 
action strength as the rapid increase in X. We did find 



that on rare occasions a further rearrangement occurs 
at larger interaction values, but it is not clear whether 
this is a genuine effect which for larger systems would 
become correspondingly less rare, or a manifestation of 
metastable configurations. 

Let us now look at the longer range behaviour of C(r). 
At half filling it has been claimed that in clean infinite 
lattice systems a second order transition to a crystalline 
state occurs at strong interactions Q . In disordered sys- 
tems, evidence of at least short range order has been seen 
in exact calculations on small systems pl|,f2l . Within the 
SCHF approximation we find no decay of correlations. It 
is well known that at half filling, nesting of the Fermi 
surface leads to a 2kp charge density wave instability, 
but away from half filling this nesting does not occur. In 
fig. 2 it can be seen that, in the presence of disorder, there 
exists no long range order in the SCHF ground state for 
ly — 1/4. This however, is also true of the related classi- 
cal system (i.e. t = and W/U ^ in the Hamiltonian 
(|l|)), where one expects the formation of a non-crystalline 
solid. One way to establish whether the electrons possess 
solid- or liquid- like correlations is to analyse the excited 
states of the system. However, to show that classical re- 
sults can provide information on the SCHF ground state 
away from half filling, at least when the particle packing 
is compact, we consider the appearance of geometrical 
frustration, where the ground state of the classical sys- 
tem contains line defects with respect to a pure crystal. 
These defects lead to the disappearance of long-ranged 
order, but at the same time give rise to magic filling fac- 
tors where (A2(A^)) exhibits large fluctuations. 

This brings us to the central observation of this study, 
namely, the strong geometry and filling factor depen- 
dent fluctuations that can arise in the average addition 
spectrum spacing (A2(A^)) as the interaction strength 
is increased [|8[. Consider first a collection of classical 
charges, on a square lattice, with nearest neighbour in- 
teractions. If the lattice is a torus with 2n x 2m sites, it 
is possible to insert up to 2nm particles without incur- 
ring an energy cost. The remaining 2mn particles cost an 
additional AU to add, so that {A2{N)) displays a peak, 
0{U), a.t N = 2nm. If on the other hand one of the sides 
(or both) of the lattice is odd (e.g. (2n-|-l)x (2m-f 1)) the 
maximum number of particles that can be added with- 
out nearest neighbours is reduced. It is not difficult to see 
that such a maximally filled configuration contains a line 
defect and long range order is lost. In other words, the 
lattice of sites without nearest neighbours is incommen- 
surate with the underlying lattice. In the commensurate 
case, the quantum system also shows a peak in A2(A^) 
at half filling, but the nesting of the Fermi surface in the 
non-interacting system also makes this filling special. We 
show below that in the incommensurate case we find that 
within the SCHF approximation, remnants of the peaks 
at the m,agic filling factors in the related classical model 
are visible far from the classical limit, despite the lack of 



nesting. 

We consider a lattice of the type (2n -|- 1) x (2n -I- 1) 
as an example, the predictions for other incommensurate 
lattices are easily obtained. In the classical limit with 
nearest neighbour interactions the first n(2n + 1) parti- 
cles can be added with no interaction energy cost, the 
next 2n + 1 particles cost an additional 2U, and the rest 
cost an additional 4t/. As a result, in the classical limit, 
(A2(n(2n + 1))) (as weh as {A2{{n + l){2n + l)))) is sig- 
nificantly larger than all other values of {A2{N)) [M. In 
our calculations we include a trivial constant interaction 
term to make the results easier to read. In Fig. S, we 
plot some typical results for (A2(iV)) for a 7 x 7 lattice, 
which shows that for [/ ^ 2 remnants of this classical 
effect can be seen clearly at the predicted filling A^ = 21. 
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FIG. 3. The addition spectrum is shown to be strongly 
dependent on the number of particles in the dot, which can 
be understood from classical arguments. The lattice is 7 * 7, 
W — 2, and the results are averaged over 400 samples. 

Similar behaviour has been observed for other sample 
sizes and geometries. Although these results correspond 
to a density regime where quantum fluctuations are pre- 
dominant, the structure in (A2(A)) agrees qualitatively 
with that of the classical counterpart. One might expect 
that extending the range of the interaction will give rise 
to a more intricate classical structure, but with corre- 
spondingly smaller amplitude which is thus more easily 
washed out by quantum fluctuations. This question is 
left for a future study. In previous work |13| the results 
of m were reproduced by interacting classical charges 
in a parabolic potential. In this case the the source of 
magic numbers incorporated the existence of topologi- 
cal defects in the ground state, as well as the interplay 
with the conflning potential. It was not a priori clear 
why such a classical model should prove useful, but our 
work suggests that quantum fluctuations strong enough 
to destablise a Wigner solid may not completely wash 
out these effects. 

In summary, we show that the metallic ground state 
develops charge density modulations, controlled by the 
electron-electron interaction, at densities r^ ^ 1 depend- 
ing on disorder. The development of the CDMs with 



increasing r^. is more rapid for short-range interactions, 
presumably because of the large gradient of the interac- 
tion potential. We also show that away from half filling, 
the CDMs are associated with short range order only. 
Elsewhere U^, it has been demonstrated that the exis- 
tence of these CDMs result in unusual fluctuations of A2 
over the disorder ensemble. Finally, we demonstrate that 
topological defects in the equivalent classical system oc- 
cur in the CDM, and that they result in strong filling fac- 
tor and geometry dependent fluctuations in A2, clearly 
visible for C/ > 2. It seems clear that the ergodicity prin- 
ciple pO| fails in this case, and so disorder averaging and 
averaging over N are not equivalent. These results lend 
support to the classical analysis of Ref. ||l£] , which sug- 
ests that the behaviour seen in the experiments of Ref. 
is due to topological defects in the classical ground 
state configuration. We stress however, that bunching is 
not generated in the geometry that we consider. 
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